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Abstract. Various coefficients of the 3.5 post-Newtonian (PN) phasing formula of 
t— 5 non-spinning compact binaries moving in circular orbits is fully characterized by the 

two component masses. If two of these coefficients are independently measured, the 
masses can be estimated. Future gravitational wave observations could measure many 
(J of the 8 independent PN coefficients calculated to date. These additional measurements 

can be used to test the PN predictions of the underlying theory of gravity. Since 
5— i all of these parameters are functions of the two component masses, there is strong 

i correlation between the parameters when treated independently. Using Singular 

Value Decomposition of the Fisher information matrix, we remove this correlations 
and obtain a new set of parameters which are linear combinations of the original 
phasing coefficients. We show that the new set of parameters can be estimated with 
significantly improved accuracies which has implications for the ongoing efforts to 
implement parametrised tests of PN theory in the data analysis pipelines. 
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1. Introduction 

Testing General Relativity (GR) and alternative theories of gravity is one of the most 
important goals of directly observing gravitational wave (GW) signals (see [1, 2] for 
reviews). The Ground based GW detectors like LIGO and Virgo are undergoing 
upgrades to second generation Advanced LIGO (aLIGO) and Advanced Virgo, making 
them capable of observing 1000 times the volume of the universe they were sensitive to 
in the first generation. The science case studies of a third generation GW detector, 
namely Einstein Telescope (ET), are being undertaken [3]. In this context, it is 
important to devise strategies which enhance the capabilities of these detectors to test 
GR and alternative theories of gravity. Inspiralling compact binaries are one of the most 
promising sources of GWs. In the adiabatic inspiral regime, the gravitational waveforms 
from these sources are accurately modelled using post-Newtonian (PN) approximation 
to GR [4]. For modelled signals, matched filtering is known to be an optimal detection 
technique in Gaussian noise. Current existing detection and parameter estimation 
techniques are based on matched filtering. Thus, prior predictability of the detailed 
nature of the waveforms becomes crucial in testing GR using inspiralling compact 
binaries. 

Since gravitational waveforms within alternative theories of gravity are not easy 
to compute, one tries to parametrise the possible deviations of alternative theories of 
gravity from GR. There have been many proposals in the literature to test GR and its 
alternatives. These proposals may be classified into two. First is to test a particular 
alternative theory of gravity by bounding the value (s) of some parameters which 
quantifies the deviation from GR. The proposals to bound scalar-tensor theories [5, 6, 7, 
8, 9, 10, 11] and massive graviton theories [12, 8, 13, 14, 9, 10, 15, 16] are examples of the 
first type. The second type is to bound the values of a set of parameters which generically 
captures the deviation from GR without referring to any particular alternative theory. 
Such a test was first proposed in [17] where the authors discussed the measurability of 
the 1.5PN term in the phasing formula. This was extended to include all the PN ordre 
terms in the phasing formula in Refs [18, 19] where the authors discussed the possibility 
of accurate determination of various PN coefficients of a binary inspiral waveform from 
GW observations and further proposed a consistency test of these phasing coefficients 
in the mass plane of the binary. In Ref [20] , a similar exercise was carried out by folding 
in the amplitude corrections to the PN waveforms. Recently there have been other 
proposals in the literature that discuss generic parametrised deviations from GR and 
how they affect GW detection and parameter estimation. Refs. [21, 22, 23] discussed 
parametrised post-Einsteinian bounds on alternative theories of gravity while Ref. [24] 
discussed implementing a variant of the test proposed in [19] in the LIGO data analysis 
context using Bayesian methods. The effect of spin precession on determining the various 
PN parameters was studied in Ref. [25] . Ref. [26] discussed possible bounds on generic 
dipolar gravitational radiation from the gravitational wave observations. 

This work can be seen as an extension of the idea proposed in Ref. [18] where 
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the gravitational waveform was parametrised in terms of the phasing coefficients. In an 
ideal scenario, one would like to treat all the eight phasing parameters as independent as 
suggested in Ref. [18]. However, since all of them are functions of binary masses, there 
are strong correlations amongst them. This results in a near-singular Fisher information 
matrix whose inversion cannot be trusted. Hence it was possible to obtain reliably the 
errors on all the parameters only for a limited range of masses, within this approach. 

In the present work, we use the singular value decomposition (SVD) approach to 
circumvent the ill-conditionedness of the Fisher matrix. Using SVD, we construct new 
parameters (corresponding to the dominant singular values of the Fisher matrix) which 
are combinations of the original phasing parameters, and obtain the errors associated 
with these new parameters. We find these errors to be very small which implies that 
one will be able to extract these new parameters with very good accuracy with ground 
based GW detectors such as aLIGO and ET. For example we focus on the first three 
dominant new phase parameters. We show that these new phase parameters of GW can 
be estimated with relative accuracy of a few times 1CT 5 as emitted by a stellar mass 
BHs at 100 Mpc. Further, for the Einstein Telescope, the accuracy for the stellar mass 
and intermediate BH binary inspirals at 100 Mpc would be 10~ 6 — 10~ 5 which would 
enable one to very precisely estimate these parameters. 

This paper is organized in the following way. Sec. 2 describes the parametrised tests 
of PN theory and discusses the waveform model we use. A summary of the noise model 
and parameter estimation scheme using Fisher matrix is presented in Sec. 3. Singular 
Value Decomposition of the Fisher matrix is explained in detail in Sec. 4. Our results 
are discussed in Sec. 5. Conclusions and future directions are elaborated in Sec. 6. 



2. Parametrised Tests of PN theory 

Fourier domain gravitational waveforms from nonspinning inspiralling compact binaries 
in circular orbit, obtained using stationary phase approximation, may be written as 

h(f)=Af- 7 <*e™W (1) 
with the amplitude A and phase ^(/) given byf 

7 

*(/) = 2nft c - $ c + j + Y,\4>k + In f}f (k - 5)/3 . (2) 

fc=0 

Here t c and $ c are the fiducial epoch of merger and the phase of the signal at that epoch. 
Chirp mass M is related to the total mass M by M = M r? 3 / 5 where rj = ^^_™^ 2 (mi 
and m 2 are individual masses of the binary) is called the symmetric mass ratio which 
is 0.25 for equal mass binaries. C which appears in the equation for amplitude of the 
waveform is a constant with < C < 1 which depends on the geometry of the source 

\ Throughout the paper we work in geometrical units with G=c=l 
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and has a root-mean-square value of 2/5 when averaged over all sky locations and 
source orientations. The coefficients in the PN expansion of the Fourier phase are given 
by [27, 28, 29]: 

^ = -^-{KMf k -^a k) ifo = -JLfrM^-^W (3) 

256 77 256 77 

where a k and a k i up to 3.5PN read as [29, 20], 

(4) 
(5) 

(6) 
(7) 

3085 2 

v + ^v 2 ; (8) 

[1 + In (6 3/2 ttM)] , (9) 

(10) 

^ _ 11583231236531 _ 640 ^ 2 6848 / 15737765635 2255 2 \ 

" 6 " 4694215680 3"^ ~ ~zT + \ 3048192 + "vi^ ) V 

76055 2 127825 , 6848, 

H ?? 2 r? 3 In (64 vrM) , 

1728 ' 1296 ' 63 V ; ' 

6848 

a 6l = - — , (11) 

/ 77096675 378515 74045 2 \ , , 

a7 = H^4016- + -J5U V - "756" ^ J ' (12) 
a fc; = 0, for fc = 0,1,2,3,4,7, (13) 

where C = 0.577... is the Euler constant. It is evident from Eq. (2) that different PN 
order terms occur with different frequency dependencies. For nonspinning binaries, each 
of the PN coefficients ijik and ipki depends only on the component masses of the binary as 
given in Eq. (3). Interestingly, from Eq. (2), the non-logarithmic term corresponding to 
k — 5 can be absorbed into a redefinition of the fiducial phase $ c as it has no frequency 
dependence. Hence there are eight parameters that characterise the phasing formula 
which are {^ , tfa, ^3, ^4, fai, ^6, ?M- 

The usual parameter estimation problem addresses the estimation of individual 
masses of the binary (or particular combinations of the component masses) using 
the 3.5PN accurate phasing formula. Two of the phasing coefficients are sufficient 
to estimate the two component masses. The additional phasing coefficients improve 
this estimation [30, 31, 29, 32]. Alternatively, one can treat each one of the phasing 
coefficients as independent parameters in the parameter estimation problem. If this is 
possible, it will enable us to test the PN theory by independently measuring each PN 
coefficient and comparing with the predictions of PN theory [18]. Though this idea 
works in principle, there will be large correlations amongst the different PN phasing 
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coefficients because all of them are essentially functions of the two component masses 
of the binary. 

To circumvent this problem, the proposal in [19] uses only three out of eight 
parameters as independent in the parameter estimation problem. This is achieved by 
re-expressing all the phasing coefficients, except the one to be used as test parameter, 
in terms of two basic parameters (usually, the two lowest order phasing coefficients). 
The third parameter used as test can be any of the higher order terms. Hence, authors 
proposed many such tests corresponding to different higher order phasing coefficients 
which are used as test parameter. Needless to say this test is not as robust as the 
original test where all of the phasing coefficients were independent, but this is one of 
the plausible ways to overcome the problem arising from correlations between various 
parameters. 

This paper paves foundation for a new approach of working with the original 
proposal where all of the phase parameters are treated independently. Using SVD, 
one can get rid of the correlations among the PN phasing coefficients and obtain the 
most dominant new phase parameters as a linear combination of the original phase 
parameters. PN theory has unique predictions about the values of these new parameters 
and this can be compared with observations by directly measuring them from the data. 
Hence these parameters can be used to do tests of PN theory. 

3. Detector noise model and parameter estimation using Fisher matrix 

We model the noise to be a Gaussian and stationary, random process. We are 
interested in the estimation of parameters of a signal buried in the noise. We use the 
Fisher Information matrix formalism [33, 34] to estimate the errors associated with the 
measurement of the signal parameters (see [35, 36] for some of the initial applications 
of Fisher matrix approach in the GW data analysis context. See a recent article about 
the caveats of using Fisher matrix as an error estimator [37]. A method to go beyond 
Fisher matrix is discussed in [38, 39].) 

Let 9 a denote the 'actual values' of the parameters and 9 a = 9 a + A9 a the best- 
fit parameters (or parameter estimate) in presence of noise. Then, in the limit of large 
SNRs, errors in the estimation of parameters A.9 = {A9 a } obey a multivariate Gaussian 
probability distribution of the form 



where M is the normalization constant. The repeated indices are summed over. The 
quantity T a b appearing in Eq. (14) is the Fisher information matrix which is given by, 



where h a = dh/d9 a . Here, (|) denotes the noise weighted inner product as defined 
below. 



p(A6) = Mexp[-T bc A9 b A9 c /2], 



(14) 



Lab = (h a | h) 



(15) 



SVD in parametrised tests of PN theory 



6 



Given any two functions a and b their inner product is defined as: 

r*wi df , OB) 

Jfrain 

where S h (f) is the one sided noise power spectral density (PSD) of the detector. Hence 
the explicit expression for the Fisher matrix is given by, 

r " b - m J f , ~m~*- (17) 

where f s is the lower cut-off frequency as given by the seismic noise limit and fi so is the 
upper cut-off frequency which is the frequency at the last stable orbit beyond which PN 
approximation breaks down. The frequency at the last stable orbit can be expressed as 
a function of the total mass M of the binary as fi so = (6 3//2 7r M) . 9?(.) denotes the 
real part of the quantity in brackets. 

The covariance matrix is the inverse of the Fisher matrix as given below 

C = {{M a M b )} = T-\ (18) 

where (•) denotes an average over the probability distribution function in Eq. (14). The 
root-mean-square error a a in the estimation of the parameters 9 a is 

a a = ((A9 a ) 2 ) 1/2 = VC^. (19) 

Note that in the above expression, index a is not summed over. 

For the noise power spectral densities of aLIGO and ET we have used the analytic 
fit given in Eq. (2.1) of [20] and Eq. (2.2) and (2.3) of [20], respectively. The lower limit 
of the frequency is assumed to be given by the seismic cut off which is taken to be 20Hz 
for aLIGO and 1 Hz for ET. 



4. Singular value decomposition and Fisher information matrix 

In this section we recast the parameter estimation problem with the help of singular 
value decomposition of the Fisher matrix. We show that the truncated SVD eliminates 
the near-singular feature of the Fisher information discussed in Sec. 2. To demonstrate, 
we opt to work in the discrete frequency domain. However, the final result does not 
depend on this choice. 

Let the 9 = {9 1 } denotes the m dimensional vector corresponding to the PN phase 
parameters as defined in Eq. (3). We divide the frequency band (f s , fi so ) into N discrete 
frequency points fk with bin-size Af such that /o = f s and /at-i = fiso- In the discrete 
domain, Eq. (17) is recast as 

Tij = 43? 

where at any frequency bin fk, h t i(fk) = dh(fk)/d9 l denotes the i th component of the 
tangent vector of the GW signal with respect to the phase parameters. 



N-l 7 



E 

fc=0 



A/. 



(20) 
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We define H = {h ik = hAM^lALx to be a matrix of size m x N where the row 

index i corresponds to the phase parameter space and the column index k is for the 
frequency index. The matrix H physically signifies the variation of the signal in the 
m dimensional parameter space 9 when evaluated at different frequency locations. It 
further acts as a building block for Fisher matrix as shown below. 
The Fisher matrix re-expressed in terms of H takes the form 

N-l 

T ij = ^2l H ikH* k + H* k H jk ] . (21) 

The main contribution to the variation of the signal in the parameter space comes 
from the variation of the phase. Hence, owing to the form of the signal expressed in Eq. 
(1), we note that H ik H* k is a real quantity in the phase parameter space 9. In terms of 
H, the m x m Fisher matrix compactifies to 

r = 2 HH t . (22) 

Henceforth, we focus on H. The singular value decomposition (SVD) is a powerful 
tool to identify the subspace of m x N space where H resides or in other words, the 
directions along which the signal varies dominantly in the parameter space at a given 
frequency location. 

The method of SVD provides a factorization of any rectangular matrix into a square 
matrix, a rectangular diagonal matrix and another square matrix as given below 

H = U H n H V^, (23) 

where Uh and Vh are the left and right unitary square matrices of m and N dimensions 
respectively and the rectangular m x N diagonal matrix fin arranges the singular values 
of H in the descending order. The columns of Uh and Vh form an orthonormal basis 
set in m and N dimensional space respectively i.e. U^U = V^V = I. Further, they act 
as (orthogonal) eigen- vectors of (Htf) mxm and (H^H)jvxat respectively. Substituting 
Eq. (23) in Eq. (22), Eq. (22) emerges as a SVD of T as 

r = u H su H ru H = su H , (24) 

where S = 2 Hh^h with the diagonal entries as = 2(f2#)| fc . 

We note that the Uh provides the SVD for the Fisher matrix with its columns act 
as eigen-vectors of the T . In other words, they provide the principal components of 
r. Hence we have shown an alternative way to obtain the principal components of 
Fisher matrix using the SVD of H. We can compare and contrast our method with the 
principal component analysis (PCA) used, for example, in Ref. [40]. 

In Ref. [40], the authors focussed on the PCA of the covariance matrix in the context 
of parameter estimation of spin modulated chirp signals in the LISA band §. Since the 
covariance matrix is simply the inverse of the Fisher matrix, both the methods would 
be equivalent for a non-singular T. 

§ Laser Interferometer Space Antenna (LISA) is a proposed space-based GW detector sensitive to the 
low frequency band between ~ 10 -4 — 1 Hz 
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4-1. Treating the ill-conditionedness of the Fisher Matrix 

To recall, in order to test PN theory, the ideal scenario would be to measure 
all the PN phase coefficients independently which was not possible due to large 
correlations among the parameters. The main hurdle to address the full problem 
was that the Fisher matrix when expressed in the 8-dimensional parameter space 
6 = {ipo,i J 2,ip3,i J i,'4'5i,4'6i,i J 6,'ip7}, is ill-conditioned. An ill-conditioned Fisher Matrix 
implies that T cannot be inverted reliably. For cases where T is near-singular, SVD 
approach provides an elegant solution to deal with the singularity and condition the 
Fisher matrix. Rest of the section, we shall focus on this. 

Consider the parameter space to be 6 = {ipo,ip2,ip3,'4'4,i J 5i,i J 6i,i J 6,i J 7} i-e- m = 8. 
We compute the SVD of H in this parameter space. In case of the near-singular matrix, 
some of the diagonal entries of £ are very close to zero. As a result, the errors in those 
parameters are very large or unreliable. Furthermore, the effect of singular nature of T 
introduces large errors in other parameters as well. Thus, to obtain the reliable errors, 
the first step is to condition the Fisher matrix or treat its singular nature. There are 
various approaches to address this problem. We follow the Truncated SVD approach in 
this paper [41]. 

The truncated SVD approach is elaborated in our context as below. We first 
perform the SVD of H. We then examine the singular values given by the diagonal 
entries of fin i.e. (flii)kk- The total number of non-zero diagonal entries gives the 
rank of H. We define the ratio of the singular values (£|ff)fcfc/(tyff)ii and consider only 
those which are greater than or equal to a small number e. Let the total number which 
satisfy this criterion be r. Then we truncate Qh to the dimensions equal to r. The 
new truncated singular matrix fij^ is r x r with diagonal entries. Consequently, the 
truncated version of left and right singular matrices and Vh is a subset of Uh 
and Vh respectively which corresponds to the first r dominant singular values. The 
truncated SVD takes the form 

H = Ukn* H v2 (25) 

where the is 8 x r, f2* is r x r and V"£[ is r x N. Thus, in truncated SVD approach 
8x N dimensional H is reconstructed by eliminating all non- dominant /near-zero singular 
values which otherwise do not contribute to H. This explicit elimination of zeros, brings 
us into a subspace in which H resides. Hence, by construction, the truncated SVD given 
in Eq. (24) provides conditioned invertible Fisher matrix 

r = u^ i s t u^ ugru^s* (26) 

where S* = 2fJ^H and u h 

give Eigen-vectors corresponding to the dominant 
Eigenvalues. Thus the truncated SVD conditions the I\ In fact, provides the unitary 
transformation Eq. (26) which diagonalises T with the diagonal entries = 2{Vl t H )\ k 
for k — 1, . . . , r. 

Comparing this approach with that of Ref. [40], one observes that the latter is 
based on the PCA of the covariance matrix which is the inverse of Fisher matrix. This 
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Figure 1. Signal to noise ratio (SNR) as function of total mass. The figure is for equal 
mass binaries. Sources are assumed to be at a fixed distance of 100 Mpc irrespective of 
the mass. A seismic cut-off frequency of 20Hz is assumed for aLIGO and 1 Hz for ET. 
Einstein Telescope has about ten times the signal to noise ratio compared to aLIGO 
for the same total mass. 



method cannot be applied if the inversion of the Fisher matrix itself is unreliable, as is 
in our case. In other words, PCA is not aimed at removing the ill-conditionedness of 
the Fisher matrix whereas our procedure is specifically devised for that. 



4-2. New Phase Parameters and the relative errors 

Let A9 a be the error in the phase parameters 9 a which obey the multi-dimensional 
Gaussian distribution following Eq. (14) as given below 

p(A0) =Jv r exp[-r a/3 A0 a A^/2], (27) 

The repeated indices are summed over and a, /3 — 1, ... 8. 

Consider the scalar T a pA6 a A6^ in Eq. (27). The truncated SVD approach should 
keep this scalar invariant as we are not removing any information but removing the 
zeros in the problem. We use this fact to obtain the new phase coordinates. 

A9 T TA9 = A6 T U^Ug A9 

= A6' T E* AO' (28) 



where superscript T shows the transpose of the matrix. The new phase coordinates 9' 
are obtained by linearly combining the PN phase parameters 9 as 

9' = Ug0 . (29) 

where 9' = {ipi,^,^, ■ ■ ■ , V'r} is a r < 8 dimensional phase vector. The Fisher matrix 
and hence the covariance matrix is diagonal in this new phase coordinates. Thus, 9' 
phase coordinates are statistically independent. 
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Figure 2. \ipk\ (left panel) and (right panel) as functions of the total mass of the 
binary for the equal mass case. 



Let the errors in parameter estimation be AO'. Owing to the statistical 
independence, the multi-dimensional Gaussian distribution of r-independent random 
variables A9' a is the product of the k Gaussian probability distributions with zero 
mean and variance o~ a = (£^ a ) -1 / 2 as shown below 

p(A9') = A/-'II Q exp[-(An7(2^)], (30) 
where a = 1, . . . r and M' is a normalisation constant. 



4-3. Phasing formula in terms of the new phase parameters 
The PN phase given in Eq. (3) can be expressed in terms of as 



8 

7T 



*(/) = 2tt ft c -®c + j + Yl P o°r ( 31 ) 

where £* =1 Pfo = El= bPk + ln/]/( fc " 5 )/ 3 such that 

P - {/- 5/3 , r\ /" 2/3 , /" 1/3 , In /, f 1/3 In /, J 1 / 3 , / 2 / 3 } (32) 

is a row vector as function of frequency. Using the unitary transformation, 9' = Ujj0 
and using UTJt = I, one can write the sum Y^j=i as 

P9 = PUU^ = P'9' . (33) 

Thus, in terms of the new phase parameters the PN phase becomes 

8 

¥(/) = 27r/f e -<Ij e + £ + (34) 

fc=i 

In the next section, we give comparison of how various behave as function of 
the total mass (for equal mass systems) in the original and new representations; see e.g. 
Fig. 2. 
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Figure 3. Left panel: Ratio of singular values of the Fisher matrix w.r.t the dominant 
one as a function of k, where k denotes the different singular values and runs 

from 1 to 8. Right Panel: For the case of aLIGO detector, variation of relative 
weights as a function of total mass for various k values. Both the panels are for equal 
mass systems (77 = 0.25). 



5. Results and Discussions 

In this section, we demonstrate the outcome of using truncated SVD approach to the 
parameter estimation problem we were interested in and present the errors on the new 
phase parameters. Before we start discussing our results it is worth recalling the typical 
SNRs expected from the aLIGO and ET detections as a function of the total mass. 
This is presented in Fig. 1. Roughly speaking, ET events will have a SNR which is an 
order of magnitude higher than that of aLIGO; thanks to the smaller seismic cut-off 
and improved sensitivity. In the figure, the distance to the source is assumed to be 100 
Mpc and we have considered binaries with equal masses. 

We now proceed and compute the Fisher matrix for a particular binary system with 
a fixed mass following the prescription outlined in Sec. 3. This amounts to computing 
the derivatives w.r.t the various parameters (the phasing coefficients) and performing 
the integration in Eq. (17). We then obtain the matrix H and carry out the SVD using 
the in-built matlab function. We observe that if we set e ~ 10~ 3 , the first 3 singular 
values are non-zero with the rest of the diagonal entries decreasing rapidly for all total 
masses considered. We retain only the first three (k = 1,2,3) leading singular values 
of S and truncate S to a 3 x 3 matrix. As was shown earlier, the above procedure is 
equivalent to first diagonalise the Fisher matrix and then rank the diagonal entries in 
the descending order and retain only the first three eigenvalues, ignoring the rest. Both 
the procedures are numerically tested to be equivalent. Values of the diagonal entries 
(normalized to 1) are plotted as a function of the index of the singular values in Fig. 3 
for two different total masses and for aLIGO and ET configurations. This figure justifies 
our omission of all eigenvalues higher than the first three. By diagonalising the Fisher 
matrix, we obtain a set of new parameters as linear combinations of the old parameters 
which are the new phasing coefficients. 
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Figure 4. Errors in the three new parameters (which are linear combinations of the 
old parameters) obtained from SVD in the aLIGO and ET bands as functions of the 
total mass of the binary for equal mass systems r\ — 0.25. Distance to the source is 
assumed to be 100 Mpc irrespective of the total mass. Einstein Telescope errors are 
better by a factor of ~ 10 due to the higher SNR ET sources have compared to aLIGO. 



In Fig. 2, we plot the old and the new |?/>fc| with the total mass of the binary. For 
the new phasing coefficients, we have shown the curves only for the dominant three 
given by the truncated SVD. The original phasing coefficients have essentially power 
law dependencies on the total mass and hence are linear in the log-log scale. The only 
exception to this is ipe which has a term proportional to M 1 / 3 In M and hence not linear. 
The new parameters which are linear combinations of the old ones are not linear in the 
log-log scale. The other interesting thing to note is that the values of the old \ipk\ 
decrease with increasing k whereas the new \ipk\ values increase with h. 

Inverse of the diagonalised 3x3 Fisher matrix S* directly gives the covariance 
matrix corresponding to the new set of phase parameters P'. The square root of the 
diagonal entries of the covariance matrix yields the 1 — a error bars associated with the 
new set of three parameters. 

In Fig. 4 we show the relative errors associated with the measurement of the three 
new parameters, as a function of mass for aLIGO and ET configurations. For masses 
< 5OM , the relative errors in the new parameters (termed as ip'^ip^ip'?) are < 10 -3 for 
aLIGO and < 10 -5 for ET. We note that removal of singularity improves the relative 
errors immensely. This is the primary result of this paper. This new set of parameters, 
which can be measured more accurately, should be more useful for devising parametrised 
tests of PN theory. 

At this stage it is worth comparing the relative errors we have obtained for the 
new phasing parameters with that of Ref [19]. In Ref. [19], the authors used a reduced 
parameter set containing two basic parameters (which are the Newtonian and 1PN 
phasing coefficients) and a third test parameter which can be any one of the higher order 
phasing coefficient. They estimated the errors associated with the measurement of these 
three parameters. The best case scenario for the third generation EGO configuration 



SVD in parametrised tests of PN theory 13 



Relative Error: Relative Error: i]/ 2 




10 20 30 40 50 60 70 80 
Total mass M (in units of Msun) 

Figure 5. Contour plots of Relative Errors in the three new parameters (which are 
linear combinations of the old parameters) obtained from SVD in the aLIGO band. 
Distance to the source is assumed to be 100 Mpc irrespective of the total mass. The 
relative errors for unequal mass cases are better than the equal mass case for fixed 
total mass. 



(somewhat similar to the ET configuration we have considered here) has relative errors 
which ranges between 1CT 2 — 1 for most of the coefficients, the only exemption being that 
of ^3 which has a relative error of 10 -3 (see left panel of Fig. 2 of Ref. [19]). Comparing 
these results with right panel of Fig. 3 of this paper, one sees that all the three new 
phasing coefficients can be estimated with relative accuracies ~ 10~ 5 which is roughly 
2-3 orders of magnitude improvement. This is solely due to the usage of the new set of 
parameters obtained using SVD. 

In Fig. (5), we plot the contours of relative errors for a general case of unequal mass 
binaries in the M — rj plane. The colors of the contours show the values of the relative 
errors. From the contours, it is clear that the approach works well for the entire range 
of rj values. Throughout the parameter space, SVD approach can provide the three 
dominant parameters (for e ~ 1CT 3 ) which can be very well estimated using matched 
filtering. To understand the relative error trends, we focus the reader's attention that 
there are two competing factors that affect the relative errors; namely the SNR and the 
number of cycles spent in the detector band. For inspiral phase, SNR oc M 5 / 6 r// 2 . The 
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Figure 6. Fixed total mass (2 M , 10 M Q , 18 M ) and varying r\ for aLIGO: (a) 
Number of cycles as a function of r\. (b) Relative error of ip[ as a function of r\. 
The analytical fit for M=2 M , error = 7.5 * lCrV' 4925 , for M=10 M , error = 
17.25 * 10~ 5 77 - 4947 and M=18 M , error = 22.8 * 1 cr 5 ?? ' 5066 - 



number of cycles N cydes oc M -5 / 3 ?^ 1 . 

For a fixed ?7 case, parameters of the binaries with low masses spend many more 
cycles in the detector band compared to the massive system. Though the SNR increases 
with the total mass, the effect of these two factors together makes the relative error 
increase with total mass as M 5 / 6 . On the contrary, for a fixed total mass and varying r\ 
case, systems with smaller r] (more asymmetric) spend more time in the band (see left 
panel of Fig. 6). Here, the number of cycles are plotted as a function of rj which varies 
as [29]. The SNR varies as 77 1 / 2 . The net product gives the relative error variation 
as r/ 1 / 2 which is inverse of the square root of the number of cycles. In the right panel 
of Fig. (6), we plot the relative error of ip[ as a function of symmetric mass ratio rj for 
fixed total mass M for aLIGO noise. Analytical fit of the relative error as a function of 
7] shows that the relative error indeed varies as ~ rf' b . 

6. Conclusions and future directions 

We revisited the parameter estimation problem where each one of the PN phasing 
coefficients are treated as independent parameters and used to devise a self-consistency 
test of PN approximation to GR. We showed that using Singular Value Decomposition 
of the Fisher matrix, it is possible to deduce the three new combinations of the phasing 
coefficients which are best estimated using matched filtering. The Fisher matrix and 
hence the covariance matrix are diagonal in this new phase parameters which in turn are 
linear combinations of the original PN phasing coefficients. Pertaining to the removal 
of the singular nature, the new phase parameters can be estimated very accurately with 
very small relative errors ~ 1CT 4 — 1CT 6 for aLIGO and ET configurations. This covers 
the BNS and NS-stellar mass BH cases for aLIGO and BNS, NS-BH and IMBH binaries 
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for ET. We studied how the errors in the new parameters vary as a function of the total 
mass of the binary as well as the symmetric mass ratio 77. 

However, here we have not discussed in detail how to rephrase the parametrised 
tests of PN theory in terms of these new parameters. This requires devising a consistency 
test using this new set of phasing parameters. One approach could be the traditional 
approach to follow the procedure outlined in Refs. [18, 19, 20] where the consistency 
test(s) where carried out in the mass plane of the binary. One may carry out similar tests 
in the mi — m 2 plane in terms of the three dominant parameters obtained using SVD. 
Another approach could be to fold in the new parameters in the Bayesian framework 
following the the work of [24] in which case the formulation and the implementation 
of the test is entirely different from the first one. Hence, our work should be seen as a 
first step towards the final goal of devising tests of PN theory where we have parameters 
which can be estimated much more accurately than those in the original proposal. These 
issues will be addressed in detail in future work. We foresee that since these new phase 
parameters pertains to the principle components of the Fisher matrix which in turn is 
related to the maximal variation of the signal, they could also be useful in the detection 
scenario. At present, we are investigating these issues. 

Note added - After this paper was finalised and circulated as a preprint, a related 
work by Brown et al [42] appeared where they addressed detection problem of inspiralling 
spinning binaries. They discuss the construction of a template bank for spinning searches 
using a waveform where all of the PN terms in the phasing are treated as independent. 
They diagonalise the metric in the 8-dimentional parameter space to obtain a reduced 
two dimensional parameter space which they use for the template placement and the 
signal detection problem. They write down the PN phasing formula in terms of a 
dimensionless variable x = f / fo where fo is some arbitrary frequency scale. The new 
parameters we defined should correspond to fo — 1 Hz in their approach. 
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